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Abstract 

Material systems with Dirac electrons on a bipartite planar lattice and possessing 
superconducting and excitonic interactions are investigated both in the half-filling 
and doped regimes at zero temperature. Excitonic pairing is the analog of chiral 
symmetry breaking of relativistic fermion theories and produces an insulating gap 
in the electronic spectrum. Condensed matter systems with such competing inter- 
actions display phenomena that are analogous to the onset of the chiral condensate 
and of color superconductivity in dense quark matter. Evaluation of the free-energy 
(effective potential) allows us to map the phases of the system for different values 
of the couplings of each interaction. At half-filling, we show that Cooper pairs and 
excitons can coexist if the superconducting and excitonic interactions strengths are 
equal and above a quantum critical point, which is evaluated. If one of the inter- 
actions is stronger than the other, then only the corresponding order parameter 
is non-vanishing and we do not have coexistence. For nonzero values of chemical 
potential, the phase diagram for each interaction is obtained independently. Taking 
into account only the excitonic interaction, a critical chemical potential, as a func- 
tion of the interaction strength, is obtained. If only the superconducting interaction 
is considered, the superconducting gap displays a characteristic dome as charge 
carriers are doped into the system and our results qualitatively reproduce the su- 
perconducting phase diagram of several compounds, like 122 pnictides and cuprate 
superconductors. We also analyze the possibility of coexistence between Cooper 
pairs and excitons and we show that, even if the excitonic interaction strength is 
greater than the superconducting interaction, as the chemical potential increases, 
superconductivity tends to suppress the excitonic order parameter. 
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1 Introduction 



In the present paper, we investigate systems of Dirac electrons on a bipartite 
planar lattice possessing competing superconducting and excitonic interac- 
tions. Excitonic pairing is the analog of chiral symmetry breaking of relativis- 
tic fermion theories and produces an insulating gap in the electronic spectrum. 
Superconducting pairing, as we know, leads to zero resistance charge trans- 
port. Condensed matter systems with such competing interactions display 
phenomena that are analogous to the onset of the chiral condensate and of 
color superconductivity in dense quark matter, which is an active subject of 
investigation pQ. The phase transitions obtained in such systems can be com- 
pared with the ones occurring in analog condensed matter systems with Dirac 
fermions. 

Indeed, several condensed matter systems have been discovered in recent years 
whose active electrons have their kinematics governed by the Dirac equation, 
rather than by the Schrodinger equation. Their dispersion relation, accord- 
ingly, has the same form as that of a relativistic particle. The reason for this 
unusual behavior of electrons whose speed is at least two orders of magnitude 
less than the speed of light can be ascribed to a particular influence of the 
lattice background on the electronic properties. 

One early example of such Dirac electrons in condensed matter can be found 
in polyacetylene [2"f5] . The undoped system is in the half-filling regime. If 
we expand the one- dimensional lattice tight-binding energy about the two 
Fermi points we obtain the linear dispersion relation associated to a massless 
relativistic particle, the two components of the Dirac field corresponding to left 
and right moving electrons. Within a large temperature range, any electronic 
excitations will have such a dispersion relation and consequently will obey the 
Dirac equation. 

In the high-Tc cuprate superconductors, Dirac points appear in the inter- 
section of the nodes of the ci-wave superconducting gap and the 2D-Fermi 
surface. Again, the low-energy excitations will correspond exclusively to these 
points jl]. 

Quite recently, it has been experimentally found that the newest high-Tc su- 
perconductors, namely the iron pnictides [5.6J also present electronic excita- 
tions whose properties are governed by the Dirac equation. Theoretical results 
also support the existence of Dirac electrons in the pnictides [7][5] . 

Another interesting example is that of graphene [9J, In this case the tight- 
binding energy corresponding to an electron in the honeycomb lattice presents 

Email address: lizardonunesOuf sj .edu.br (Lizardo H. C. M. Nunes). 



2 



a band structure such that the valence and conduction bands touch precisely 
in the vertices of two inequivalent Dirac cones and any electronic excitations 
appearing in the conduction band will have the dispersion relation of a rela- 
tivistic massless particle and their properties, accordingly, will be determined 
by the Dirac equation. In the half-filling regime, with a completely filled va- 
lence band, we will have just a Fermi point at the vertex, whereas in the 
presence of doping, a Fermi surface builds up. 

For many technological applications, it would be interesting if graphene were 
a semiconductor, instead of a semimetal. Therefore, the possibility of an insu- 
lator state or the appearance of an excitonic gap due to short range electronic 
interactions in a system of massless Dirac fermions on the honeycomb lattice 
has been the object of intense investigation [TU|11|I12|I13|I14] . Nonetheless, early 
proposals that electron-electron interactions could generate an electronic gap 
were investigated even before the discovery of graphene [ToTTE] , where the gap 
opening is an analogue of the "chiral symmetry" breaking process that occurs 
in massless quantum electrodynamics (QED) in two dimensions [?{H]. This 
is the type of excitonic instability that we consider in the present paper. 

The fact that some of the above materials become superconductors upon dop- 
ing provides a vast phenomenological motivation for the investigation of the 
superconducting phase diagram of Dirac electrons systems. Therefore, we anal- 
yse the effect of Cooper pair formation in such systems, paying special atten- 
tion to the cases where doping is present, in which case we must introduce a 
finite chemical potential. 

Indeed, we consider a quasi-two dimensional model consisting of Dirac elec- 
trons both with superconducting and excitonic interactions. In particular, we 
study the competition between these two interactions and the resulting differ- 
ent situations that arise as we independently vary the coupling parameter of 
each interaction term and as we add charge carriers to the system. 

The paper is organized as follows: In Sections [2] and [3] we obtain the continuum 
limit of our model and we calculate the free-energy (effective potential) at zero 
temperature. In Section H] we analyze the phase diagram of the model in the 
absence of doping or impurities, in Section [5]the effect of doping is described by 
the introduction of a non- vanishing chemical potential and the physical effects 
of the interplay between the amount of doping and the interactions strength 
bring interesting consequences, which we fully analyze. Each interaction is 
firstly taken independently and its respective phase diagram is obtained. At 
the end of Section [5] we combine our results and we investigate the possibility 
of coexistence between Cooper pairs and excitons. 
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2 The model 



Consider a system presenting a bipartite lattice formed of sublattices A and 
B, such as in graphene [9] or in the pnictides [6]. Let a] CT = J2k e * k ri a k,o- an d 
fc| CT = J2k e * k ri o- De ! respectively, electron creation operators, with spin a, 
on site % of sublattices A and S. We take the non- interacting hamiltonian as 

m 

k,cr 

-t £ s k (a k ,A,<x + h.c. ) , (1) 

k,<r 

where the /i is the chemical potential and the last line in the RHS of the above 
equation describes the hopping of electrons between adjacent sites of different 
sublattices A and B. 

For the case of graphene, we would have t ~ 2.8 eV and, for the corre- 
sponding honeycomb lattice, = 1 + e ikai + e* ka2 where &i = ae x and 
2a2 = a (e x — V3e y ^j , with a as the lattice parameter. 

We assume the interaction has two terms leading, respectively, to supercon- 
ducting and excitonic instabilities. The isotropic s-wave pairing state is de- 
scribed by a BCS-type interaction as 

-f^sc = ~Kc ( a k,(7 a ~k,-CT a -k',-o- a k',cr 
k,k',cr 

+ OV^-^v) > ( 2 ) 

The value of A sc is to be determined by some underlying microscopic theory, 
which is not considered here. In the present paper, the strength of the super- 
conducting interaction is a free parameter of the hamiltonian and can take 
any arbitrary positive value. 

The superconducting order parameter, by definition, is given by 

- A = A sc £<4,/-k,;> = A sc X^-Ja ) • ( 3 ) 

k k 

Using a mean-field approximation, we insert the above expression in Eq. Eq. (j2J), 
thereby rewriting the superconducting interaction as 

#r = e a Ry.k,-. + 0- k ,-.) + h - c - > ( 4 ) 

k 
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except for a constant term. 



We define the four-component Dirac spinor as ^/ r a = (a T>(r , b r a ) [9]. In the 
spinor representation, the Coulomb interaction for an unitary volume reads 

ff C = EE *r,„ 7° ®r,a U(t, t') ^ 7 ° %>,*> , (5) 
cr, a' r,r' 

where \j/ r(T = \I/J CT 7 , with the Pauli matrix 7 , and U(r,r') is the Coulomb 
potential. A sufficiently strong Coulomb repulsion may set the stage for an 
excitonic instability which results in the opening of a gap in the electronic 
spectrum. The transition is signaled by the development of the excitonic order 
parameter [15] 

a = = A cxc E(«L«k, CT - &L,A,*> , (6) 

which is identified as a formal analog of the chiral symmetry breaking in 
relativistic fermion theories. This is the type of excitonic instability that we 
focus upon below. In momentum space, the corresponding excitonic pairing 
interaction in a mean-field approximation becomes, 

=£*(4^-^A*) • (7) 

k,cr 



Combining it with Eq. Eq. and Eq. Eq. (j4| 
tonian, except for constant terms, 

k 

where the matrix Ak is given by 



we obtain our model Hamil- 



(8) 



Ax 



[i + cr —tSk 
-ts* k /i -a 
A* 




A 

a — ji 



A 


tsl 



V 



A* 







ts k -(11 -a) j 



and 



(9) 



In the case of graphene, there are six Dirac points at the corners of the first 
Brillouin zone in the honeycomb lattice; however, only two of them are non- 
equivalent. We denote them by K = — 4ir/3ae x and K' = 47r/3ae x and we 
expand Eq. Eq. (jSJ) in their vicinity. Collecting only the linear terms for Sf., 
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and taking the continuum limit of our model, we obtain 

= E / 7^2 Aa *„(*) - ^ - ^ , (10) 

a (Z7TJ A sc A exc 

where a = K, K' denotes the Dirac points and the matrix A a is given by 
Eq. (jHJ), replacing tsk by — vf (fez — ik y ) or — t>p (A»e + ik y ) in the vicinity of K 
or X'. Moreover, the Fermi velocity is given by Vp = y/3ta/2, with h = 1 from 
now on. 



3 The effective potential 



The partition function in the complex time representation can be written as 



Z = — I £>II£>$ 

Zq 



x exp 



where $ is an arbitrary set of fields, II is their conjugate momentum and Zq, 
the vacuum functional. 

Inserting Eq. ( TTOl) in the above expression, we can perform the gaussian in- 
tegral over fermionic fields to obtain an effective potential. Conversely, by 
integrating on A and a we would reobtain the original fermionic quartic in- 
teraction hamiltonian by means of a Hubbard-Stratonovitch transformation. 



The effective potential per Dirac point reads 



IAI 2 a 2 



cff : 



CO 

-t E 



sc /v exc 

00 ' d 2 k 



x log 



[J 



m=i [ioj n -E j (a = 0,A 



(12) 



where o; n = (2n + l)7rT is the Matsubara frequency for fermions (with the 
Boltzmann constant ks = 1) and the four Ej's are 



S.- = ±W|A| 2 + 



\Jcr 2 + (v F k) 2 ± /i 



(13) 
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Notice that Eq. (|T3|) is the typical dispersion relation observed in effective 
models for quantum chromodynamics, where E~ corresponds to the energy 
required to create a particle or a hole above the Fermi surface and Ej~ is 
the corresponding antiparticle term [2D]. Moreover, the effective potential in 
Eq. (Tl2]) is the same obtained in [21], where the formation of chiral and di- 
quark condensates involving two quark flavours in QCD was investigated. In 
their case, the values for the quark-antiquark and quark-quark interactions 
strengths were set equal by construction, via a Fierz transformation applied 
to the Nambu-Jona-Lasinio model, which is not the case in the present paper, 
since A sc and A exc are taken to be free parameters. 

At zero temperature, Eq. (fT2l becomes 



cff : 



|A| 2 

A sc 
d 2 k 

W~ 2 



+ 



a 

Aexc 

00 dw 



2tt 



log 



7 (lo, k, A, a) 



7(w, k, A = 0,ct = 0) 



(14) 



where 



7 (lo, k, A, a) = (| A| 2 + a 2 + v\k 2 + u 2 f + /i 4 
+2(\A\ 2 -a 2 -v 2 F k 2 + cu 2 ) fi 2 . 



(15) 



Calculating the integral over u in Eq. ( |T^|) . the effective potential becomes, 



IAI 



a- 



eff ■ 



A sc A c 

d 2 k 



J2J\A\ 2 + tf-2v F k 



L / 



(16) 



where ^i(x) = \/x + a 2 + Z/x, with I = +1 or —1. Introducing the momentum 
cutoff A/vp, we arrive at the following expression for the effective potential, 



_ j. |A| 2 a 2 
VcS= A + A~ 



2A^ 



-^E[£?(a 2 )-£?(o) 



3a 3a 



+ 2^£'M{6(*)^(s) 



+iA| 2 io g [e/(x) + ^(^)]} 



(17) 



where a = 27rt>p and -Ez(x) 



A| 2 + ef(x). 
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In particular, for // = a = 0, Eq. ( ITT)) is reduced to 



off 



|A| 



2 

3a 



A| 



A 



3 

2\ 2 



IAI 



A 1 



which is exactly the same result obtained by some of us for a theory describing 
a single layer of Dirac electrons interacting via a BCS-type superconducting 
interaction. Analysing the minima conditions for the above effective potential, 
we find a quantum phase transition at the critical coupling A c = a/A [19J , 
leading to the onset of superconductivity in the system. 



In the remaining, we analyze the conditions for the appearance of supercon- 
ductivity or excitonic fluctuations at zero temperature, taking into account 
the two competing interactions present in our model. 



4 Zero chemical potential 



We analyze in this section the zero temperature phase diagram of our model 
in the absence of doping or impurities, namely, at zero chemical potential. We 
search for values of A and a that minimize the effective potential. The con- 
ditions for the appearance of superconductivity or the excitonic fluctuations 
are furnished by the existence of nonzero solutions of the respective order 
parameters, which minimize the effective potential. 

The results are obtained in terms of an arbitrary energy cutoff A (momentum 
cutoff A/vp), which is always provided by the lattice in condensed matter sys- 
tems. Indeed, we have A = 2nhv-p/a. Since a is the smallest distance scale, A 
becomes a natural high-energy cutoff. This frequently happens in condensed 
matter. A familiar example in the case of conventional, phonon mediated su- 
perconductivity, is the Debye frequency (energy) a natural cutoff that emerges 
in BCS theory. 

For \i = 0, Eq. (fTTj) becomes 



cff : 



+ 



2 

3a 
IAI 2 



A| 2 + a 2 + A 

.2 



3 

2\ 2 



A| 



a 



2\ 2 



A^ 



+ 



<7 



A S r A, 



(19) 



which can be rewritten as 



8 



(A 



(20) 



where a = oj A 2 and A = |A|/A 2 . 

2.0 r 




Figure 1. Order parameter £q as a function of the interaction strength 
A — A sc — A eX c- 



First, let us assume the following situation: A = A sc = A cxc . In such case, we 
can define a new order parameter, ( 2 = |A| 2 + cr 2 , which is inserted in Eq. (119j) . 
and we get exactly the same effective potential given by Eq. ( 1T81) . Therefore, 
analyzing the values of £ that minimize the effective potential, we arrive at 
the following, 



A | 2 + ^ 



, A < A c 



(21) 



where £ , A and a denote the solutions for the minimum of V e g. The plot of 
the order parameter as a function of the interaction strength can be seen in 
Figure [TJ As in the case when we had just a superconducting interaction and 
H = 0, we find a quantum phase transition at the critical coupling A c = a/A. 
However, we now have also the possibility of coexistence between supercon- 
ductivity and the excitonic fluctuations if the interaction strengths are equal. 

Indeed, as can be seen in the tridimensional plot shown in Figure EJ the minima 
of the Kff presents a radial symmetry, in agreement to the result in Eq. (l2Tj) . 
Therefore, it is possible to find simultaneously nonzero values for a and A 
which minimise the effective potential, indicating that excitons and Cooper 
pairs coexist in the system. 
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0.7 0.7 



Figure 2. Effective potential including both cases: excitonic and superconducting 
interactions. The parameters are taken as fi = 0.0, A exc = 1.5 and A sc = 1.5. All 
quantities are given in units of A. 

The same result was obtained in the framework of nuclear physics [2T] , where 
the formation of chiral and diquark condensates involving two quark flavors in 
QCD was investigated. In their case, the values for the quark-antiquark and 
quark-quark interactions strengths were set equal by construction, via a Fierz 
transformation applied to the Nambu-Jona-Lasinio model. 

In our case, on the other hand, we are not restricted to the relation A sc = A exc 
and we shall analyze in what follows the situation for which A sc ^ A cxc - 

We first assume that A sc > A exc . Then, Eq. (12"0"|) can be written as 



V, 



cS 



A 3 



Ah 



(c 2 + i)"-c 3 -i 



where 



with A P 



^ _ A sc — A c 

■^sc ^exc 

A exc / A c and A sc = A sc / A c . 



> 0. 



(22) 



(23) 



Analysis of the second derivatives of the effective potential yields the minimum 
of Veg at the point 

fXl-1 



(Aq, ct ) 



2A S 







(24) 



when A sc > 1 and 5 > for positive values of Co an d ct . Hence, Co — ^o an d 
we have only superconductivity in the system provided A sc > A exc . Indeed, 
Figure [3] shows the plot of V e s for A sc > A exc . There is a single minimum for 
positive values of A and a exactly at the point given by Eq. (T5] 
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Figure 3. Effective potential including both cases: excitonic and superconducting 
interactions. The parameters are taken as \i = 0.0, A cxc = 1-5 and A sc = 1.75. All 
quantities are given in units of A. 

Accordingly, when A sc < A exc , the minimum becomes 



when A cxc > 1 and we only have excitons in the system. Figure H] shows the 
plot of V e g for A sc < A cxc and the single minimum is given by Eq. ( |25l) . 

In summary, in the undoped regime, where \i = 0, the above results show that 
we can have coexistence of excitons and Cooper pairs whenever the corre- 
sponding interactions strengths are equal. If one of the interactions is stronger 
than the other, then only the corresponding order parameter is non- vanishing. 



5 Finite chemical potential 

Consider now the effective potential in Eq. ( Fl6|) . in the case of a nonzero 
chemical potential. Making the change of variable x = (vpk) 2 , we get 




(25) 



0.7 




Figure 4. Plots of effective potential including both cases: excitonic and supercon- 
ducting interactions. These plots are the same but seen of different points of view. 
Parameters are taken as /j, = 0.0, A exc = 1-75 and A sc = 1.5. All quantities are given 
in units of A. 

|A| 2 a 2 2A 3 



V e 



cff ■ 



A sc A cxc 3 a 

I f A 2 



-— W dx \/l A l 2 + (V^T^+ln) 2 . (26) 

We notice that V e s(fi) = V e s(— /i). Therefore, given A sc and A exc , the minima 
for fi > are exactly the same for fi < and we constrain our analysis to 
positive values of the chemical potential in this section. 

In the following, we consider superconductivity and excitonic fluctuations sep- 
arately, and then we combine our results to study the superconducting and 
excitonic interactions together as the chemical potential increases. 



5- 1 A exc 7^ 0, X sc — 



In the absence of the superconducting interaction, the superconducting order 
parameter is zero and Eq. ( 1261) becomes 



12 



a 



0) 



a 2 



A, 



, Jo 



y + a 2 + I ft 



+ 



(27) 



where A exc = A exc / A c and the quantities within the modulus are divided by A 
and are denoted by a and ft, since we have made the replacement y = x/A 2 
in the above integral, as in the previous section. 



Since \i > 0, we have that y/y + a 2 + /i = yjy + a 2 + ft, for every positive 
values of a, ft, y E [0, 1] and I in the above equation. For / = —1, on the other 
hand, according to the modulus definition, we have 



y + a 2 - ft 



y/y + a 2 - /j, ,ft<Vy + & 2 



yy + a 2 - ft J , ft > y/y + a 1 
Hence, we have three different results for the effective potential: 
• ft < a: 

In this case, the modulus becomes 



(28) 



y + a 2 - ft 



y + d 2 -ft, Wye [0, 1] 



(29) 



and Eq. ([27]) yields 



V. 



(i) 



a 2 2 



eil' 



A fixr 3 



a 2 + lY -a 3 -l 



(30) 



which is exactly the same expression given by Eq. (1181) if one replaces a and 
A exc for ( and A sc respectively. 



• ft > a/1 + 5- 2 : 

The modulus becomes 

and Eq. (EJJ) yields 



Pr2 



-Uy + d 2 -ft), Wye [0,1] 



v£ = ^-ft + 



A,, 



(31) 



(32) 



• a < fi < Vl + cr 2 : 

At this intermediary condition, Eq. ( 1271) yields 
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(3) 



a 2 



eff 



A 



cxc 

1 

2 



dy [\ly + a 2 - fi 

+ f dy (v^-/i 



(33) 



and this equation becomes 



(3) 



a 2 



eff 



a 2 + l A§ 



A* ~2 



(34) 



In the first case, the function is minimized by 



Aexc < 1 
; A exc ^ 1 



(35) 



Notice that this result does not depend of a given chemical potential, since the 
effective potential is /i-independent . In the second case, V^ s is minimized by 
a = 0. In the third case, we do not find a nonzero minimum for V^ s as well, 
but only a maximum at the physical value of \J (A -1 + /i) 2 — 1. Therefore, the 
system does not present excitonic fluctuations when jl > y/1 + a 2 or a < jl < 

Combining all the above results, we can obtain the phase diagram of the 
system for the excitonic fluctuations. 

Let us start by the case A exc < 1 and jl = 0, in this case, ct = 0, according 
to Eq. ( 135]) . Therefore, as the chemical potential increases, we can only get 
the second or the third situations for the effective potential, V^ ff or V e g re- 
spectively. Hence, we never reach a nonzero minimum for the excitonic order 
parameter whenever A cxc < 1 for all values of jl. 

When Aexc > 1, the analysis is subtler, but also straightforward. We start 
by the assumption that jl = 0. In this case, we do have a nonzero minimum 
ct given by Eq. ( 133]) . Now, if we take a small value for jl, say jl = jl, as 
a increases, the effective potential passes from the second situation, V^ s , to 
the first situation, , at the point a = jl and the effective potential can 
only reach a minimum if do > jl as a continues to increase. However, even if 
ct > jl, the effective potential may or may not possess a minimum, because 
the minimum for v}^ has to be smaller than any value of V}^\ when a < jl. 

~ ~ (2) 

Since A exc > 1, we have that V^ s is a monotonically increasing function for 
all ct < ji. Therefore, we get the following condition for the appearance of a 
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minimum for the effective potential: V}q (a = <To) < V^'(<j = 0). So, whenever 
this condition is satisfied, we have a nonzero excitonic order parameter, a = 
(A 2 - l) /2A, for A exc > 1 and jl > 0. 



Moreover, notice that the particular case when v}^\a 
provides a critical value for the chemical potential, 



*„) < V$(a = 0) 



+ - 



+ 2 



XL 




(36) 



As /2 increases above the threshold ft > fi c , we can no longer find a minimum 
for the effective potential. Hence, provided A CX c > 1, ft = He is also a quantum 
critical point. 




Figure 5. Phase diagram of the system taking into account only the excitonic 
interaction. The gray area indicates the the presence of a nonzero excitonic order 
parameter in the system. The solid curve in the plot is given by the critical chemical 
potential in Eq. (i36l) 



Figure |5] shows the phase diagram as a function of A exc and jl. The area below 
the curve for /i c (A exc ) indicates the region where the system presents excitonic 
fluctuations. To our best knowledge, this is the first time that this phase 
diagram was obtained. 
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5.2 X sc 0, X exc — 



In this section we analyze the phase diagram of the system taking into account 
only the presence of superconducting interactions. Hence, in the absence of 
excitonic interactions, Eq. (126]) becomes 




-^E/ Q 1( Wl A l 2 + (V^+^) 2 , (37) 

where the relevant quantities are divided by A, as in the previous sections. This 
problem has been previously studied by some of us for // = [H], and some 
results for fi ^ have been previously reported considering the s-wave and an 
exotic p-wave pairing [22]. For /i = 0, the system undergoes a quantum phase 
transition at the critical coupling A sc = a/A = A c and the system becomes 
superconducting whenever A sc > A c . 

Given values for jl and A sc , we can look for positive numerical values of A that 
minimise the effective potential, given by Eq. (|37|) . Once again, our analysis 
is constrained to positive values of /i, since V e ^{jl) = V e s(— jl). 

0.6. 1 1 , 1 1 1 , 1 1 1 1 1 1 , , 



0.4 



0.2 



0,0 1 1 1 1 1 1 1 1 1 1 1 1 J — i 1 J 

0.0 0.5 1.0 1.5 

Figure 6. Plot of the superconducting gap |Ao| as a function of the chemical po- 
tential jl. The superconducting interaction is taken to be A sc = 1.5. 

We start our discussion considering the case A sc > 1. As can be seen in Fig- 
ure El for jl = 0, the minimum of the effective potential A is given by Eq. ( )2T|) 
implying ( — > A . As jl increases, A reaches a maximum at an optimum value 
of the chemical potential and A decreases as jl increases even further. 
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Figure 7. Plots of the superconducting gap Ao as a function of the chemical potential 
jl for several values of A sc < 1. 

For the case A sc < 1, our results are shown in Figure [0 Starting at jl = 0, the 
system is in the normal state. However, as ft increases, the system asymptoti- 
cally becomes superconducting and the order parameter also increases up to a 
maximum value at an optimal chemical potential. As jl increases even further, 
A decreases and the energy gap displays a dome-shaped plot. 

Our results are consistent with [23J , where chiral and diquark condensates are 
calculated for two-color and two-flavor QCD. For a certain choice of param- 
eters, the authors obtain the same dispersion relation given by Eq. ( fl~3l) and 
their numerical results for A also display a dome-shaped plot, as can be seen 
in Fig. 1 of [23], for the choice of parameters I, referred as the weak-coupling 
case. 

Since the energy gap and the superconducting critical temperature T c are 
proportional, our results qualitatively reproduce the superconducting phase 
diagram of several compounds, like 122 pnictides and cuprate superconductors, 
where the critical temperature displays a characteristic dome as charge carriers 
are doped into the system. 

A dome-like structure of the superconducting phase for strongly interacting 
two-dimensional Dirac fermions has been previously obtained in [21], where 
the authors observed a dome-shaped plot of the superconducting phase at 
intermediate filling fractions, surrounded by the normal phase for fillings close 
to unity or zero, which is consistent to our results. 

We conclude this section with a final remark: contrary to what happens at 
H = 0, where superconductivity appears in the system only when A sc > A c , at 
finite chemical potential, we could always find a finite jl that provided nonzero 
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superconducting gaps, even at small values of A sc , as can be seen in Figure [7J 
where plots of A as a function of p, for small values of A sc are shown. This 
result suggests that even for very small superconducting interaction strengths, 
one can always find superconductivity as charge carriers are added to the 
system. 



5.3 A sc ^0,A sc ^0 



In the previous sections we have analysed the conditions for the appearance of 
the excitonic and superconducting order parameters considering each interac- 
tion separately. In this section, we study the minima of the effective potential 
taking into account both the interactions simultaneously and also the variation 
of the chemical potential. From Eq. f )26|) . V e s can be written as 



V cS =—V cS (X sc ^0, A sc ^0) 
|A| 2 a 2 2 



-\ E f d y y I A l 2 + (Vv + ** + 2 (38) 

and, once again, we constrain ourselves to positive values of the chemical 
potential, since V e s is even with respect to /i. 

We start our analysis with p, = 0, assuming that both A exc , A sc > 1. We 
also impose that A exc is larger then A sc . Therefore, as we have seen in Sec. 
HI since the interaction strengths are different we cannot have coexistence of 
superconductivity and excitonic fluctuations. Actually, since A cxc > A sc , we 
have zero superconducting order parameter and a = (A^ xc — l)/2A exc . Indeed, 
as can be seen in Figure HI the point (0, do) is the minimum of the effective 
potential for non- negative values of do or Ao- Moreover, as we start to increase 
the values of the chemical potential, the minimum of the effective potential 
remains the same as long as we constrain ourselves to small values of ft, as 
can be seen in Figure [SJ 

Interesting results appear when jj, continues to increase, since we regain the 
radial symmetry of the minima for the effective potential, as can be seen in 
Figure [H In other words, we recover the coexistence of superconductivity and 
excitonic order fluctuations even if A cxc > A sc as the chemical potential reaches 
a certain value. 

Furthermore, as the chemical potential increases even more, other interesting 
results show up, as can be seen in in Figure [TUl In that case, the plot for the 
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-0.7 0.7 




Figure 8. Effective potential including both cases: excitonic and superconducting 
interactions. The parameters are taken as jl = 0.1, A exc = 1.75 and A sc = 1.5. All 
quantities are given in units of A. 




Figure 9. Effective potential including both cases: excitonic and superconducting 
interactions. The parameters are taken as jl = 0.3, A exc = 1-75 and A sc = 1.5. All 
quantities are given in units of A. 

effective potential shows the minimum at the point a = and finite A , 
which means that the excitonic fluctuations were suppressed from the system 
and only superconductivity remains. In other words, even if A exc > A sc , as 
charge carriers are added to the system, doping effects eliminate the excitonic 
order parameter and favours superconductivity. 

Our results are consistent to the results obtained in hadronic physics, where, 
as the baryonic chemical potential increases, chiral symmetry in the system is 
restored and color superconductivity sets in pp. 
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a 0.7 0.75 

Figure 10. Effective potential including both cases: excitonic and superconducting 
interactions. The parameters are taken as fi, = 0.5, A ex c = 1-75 and A sc = 1.5. All 
quantities are given in units of A. 

To conclude this section, we analyze the case when both A cxc , A sc < 1. At 
zero value of the chemical potential, we do not have the presence of super- 
conductivity or excitonic fluctuations in the system. However, as can be seen 
in Figure HTJ the plot of effective potential show that we can find a nonzero 
superconducting gap for, A cxc = 0.9 and A sc = 0.8 at jl = 0.5 Therefore, 
even for small interactions with A sc < A cxc , there is the possibility of finding 
superconductivity as charge carriers are added to the system. These results 
indicate that, as the chemical potential increases, the insulating gap in the 
quasi-particle dispersion goes to zero, and that superconductivity sets in. 



6 Conclusions 

Our investigation has produced several different results concerning different 
regimes of the coupling parameters of the superconducting and excitonic in- 
teraction terms on materials possessing Dirac electrons. We also investigate 
the effects of doping on these. 

At fi = 0, we have shown that Cooper pairs and excitons can coexist if the 
superconducting and excitonic interactions strengths are equal and above a 
quantum critical point, A = A sc = A cxc > a/A. If one of the interactions 
is stronger than the other, then only the corresponding order parameter is 
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Figure 11. Effective potential including both cases: excitonic and superconducting 
interactions. The parameters are taken as jl = 0.5, A CX c = 0.9 and A sc = 0.8. All 
quantities are given in units of A. 

non- vanishing and we do not have coexistence. 

For fi 7^ 0, taking into account only nonzero values for the excitonic interaction 
strength, a critical chemical potential, as a function of A exc , was obtained, 
as shown in Eq. (|36|) . Given A exc , values for the chemical potential below 
Hc(Kxc) indicate the presence of excitonic fluctuations in the system. To our 
best knowledge, this is the first time that this phase diagram was obtained. 

If there is only the superconducting interaction in the system, a remarkable 
result was obtained for the case: A sc < a/ A the superconducting gap displays 
a dome-shaped curve as a function of the chemical potential, indicating the 
appearance of superconductivity as charge carriers are added to the system. 
We emphasize that this curve actually vanishes exponentially as fi — > 0, hence 
superconductivity persists down to \i = 0, below some temperature. This is in 
agreement with the fact that Cooper's theorem must be valid for /x ^ 0, when 
a Fermi surface builds up. 

Finally, we have analyzed the possibility of coexistence between Cooper pairs 
and excitons for /i ^ and we show that, even if the excitonic interaction 
strength is greater than the superconducting interaction, as the chemical po- 
tential increases, superconductivity tends to suppress the excitonic order pa- 
rameter, even if A exc > A sc . 

We have considered here exclusively the zero temperature situation. We are 
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presently investigating finite temperature effects on this system and will report 
the results elsewhere. 
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